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We show that a single impurity embedded in a cold atom bosonic Mott insulator leads to a novel 
polaron that exhibits correlated motion with an effective mass and a linear size that nearly diverge 
at critical value of the on-site impurity-boson interaction strength. Cold atom technology can tune 
the polaron's properties and break up the composite particle into a deconfined impurity-hole and 
boson particle state at finite, controllable polaron momentum. 

PACS numbers: Valid PACS appear here 



The exploration of and unprecedented control over 
quantum many-body systems have become central driv- 
ing forces in cold-atom physics. Optical lattices (the 
standing wave patterns of frequency-stable, reflected 
laser beams that are experienced by ultra-cold atoms as 
periodic potentials have unlocked strongly-correlated 
lattice physics to cold atom simulation. The Bose- 
Hubbard Hamiltonian ^] has successfully and quanti- 
tatively modeled boson atom dynamics in such optical 
lattices. Studies of the boson superfluid to Mott insu- 
lator (MI) phase transition predicted within the Bose- 
Hubbard model Q highlight the unusual cold atom ac- 
cess and control: experiments observed the transition 
revivals of inter-site superfluid coherence [5.] and differ- 
ent integer filling-number islands in the MI phase as seen 
spectroscopically |6| and, more recently, by direct imag- 
ine: In this Letter, we show that by combining the 
control over optical lattice parameters (varying the bar- 
rier height), over the inter-particle interactions (varying 
an external, homogeneous magnetic field in a Feshbach 
resonance [il HI), and over particle-species (creating 
mixtures of distinguishable kinds of atoms), cold atom 
experiments can realize a novel [l8| and controllable po- 
laron in the MI phase. 

The polaron state is induced by an impurity atom that 
experiences the same (or similar) optical lattice potential 
as the MI bosons from which it is distinguishable. The 
polaron consists of the impurity and a boson that is pro- 
moted to the next Hubbard band by a strong impurity- 
boson repulsion. If the impurity-boson interaction is Fes- 
hbach tuned to be nearly as repulsive as the boson-boson 
interaction, the excited boson remains pinned to the im- 
purity site. The boson-impurity pair propagates through 
the lattice with finite total momentum K. The linear 
polaron size, Ak, or average distance between the impu- 
rity and the excited boson sensitively depends on K, and 
increases with increasing impurity-boson repulsive inter- 
action UiB [HI- Surprisingly, the effective mass of the 
polaron increases with its size. Moreover, in the strong 
coupling limit the effective mass diverges at the criti- 
cal value oi UiB above which the boson-impurity pair 
becomes unbound (Aq oo). Most significantly, opti- 
cal lattice experiments can create probe and manipulate 
composite particles with properties that are unusual in 



traditional polaronic and excitonic systems. 

We will model our problem with a bosonic Hubbard 
Hamiltonian on a hyper-cubic lattice of dimension d: 

T-L = —ts ^ hlh^> —tj ^ cjcj./ 

(r,r') (r,r') 

+ ^ E ^rblbrbr + UjB ^ blbrcU^. (1) 
r r 

The operator bl (6^.) creates (annihilates) a boson on site 
r, while cj (c^.) creates (annihilates) the impurity, (r, r') 
indicates that r and are nearest-neighbors. The statis- 
tics of cl and c^ is irrelevant because we are considering 
a single impurity problem. 

Here we will only consider the case of strongly repul- 
sive on-site boson-boson interaction, Ubb ^ ^|^bU^/|, 
and integer filling factor (6j6r) = n to stabilize the MI 
state. An increase in the intensity of the optical lattice 
laser of wavelength A so that the lattice height Vo signif- 
icantly exceeds the recoil energy Hujr (= /i^/pm^A^]), 
where is the actual boson mass) drives the boson 
system deep into the MI regime. This increase tightens 
the trapping frequency ujt of the optical lattice wells, 
ujt = 2ujji^, where x = Vo/huR. For fixed boson- 
boson scattering length a^B, an increase in x, x > 1, en- 
hances Ubb ^ ficjT>/87r (<^bb/A) x^/^, and exponentially 
decreases the hopping matrix element tBB which we es- 
timate as ^ {3/2)hwTe~^^^'^^ ^ . Scattering lengths 
take the value of a few nm whereas the optical wavelength 
is of order of a micron but even for x ~ 9, [/^b can ex- 
ceed tB by an order of magnitude whereas t b can still be 
of the order of 10 Hz. Thus, the Ml-regime can be ac- 
cessed by varying the optical lattice height, which leaves 
a homogeneous magnetic field tuned near a Feshbach res- 
onant value as a control knob to vary aBB (oc Ubb) or 
the impurity-boson scattering length aiB (cx Uib) inde- 
pendently. We consider the regime Uib ^ ^|^bU^/|- 
The eigenstates of 1-L become highly degenerate in the 
static limit tB = tj = 0. The lowest-energy eigenstates, 
illustrated in Figdl can be expressed as 

|R,r) = J_^ clbl^^b^\0) for r ^ 0, (2) 
Vn(n + 1) 

|R,0) = 4|0), (3) 
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where |0) = [Ir and where |0) denotes 

the vacuum state. The states |R, r) correspond to the 
impurity at site R creating a hole at the same site and 
relocating the removed boson to site r (see Fig. HJa)). 
The size of the resulting polaron is the size of the particle- 
hole bound state, i.e., the mean value of r. 

The energy eigenvalues of the states (|3j) are given by 



E - Uss+{n-l)Ujs+E, 



'H|R,r) = {UBB^UiB{n-l) 
?^|R,0) = (n/7/B+Eo)|R,0) 



•£;o)|R,r) for r^O 
(4) 



where Eq = NUn{n — l)/2 represents the ground state 
energy of the undoped MI: n\0) = Eo\0) {N is the total 
number of lattice sites). The energy of any other eigen- 
state is of order Ubb or Uib higher than the energy of 
the eigenstates (j4]). According to perturbation theory, to 
lowest order in the small parameters tj^/llBB and tj^/UjB 
{v = I^B)^ we can exclude those high-energy states and 
restrict the action of % to the lowest energy subspace 
generated by the states (|4]). In this way we reduce the 
many-body problem to a two-body system described by 
the low-energy effective Hamiltonian, 



n = ^HoK + HlK, 
K 

noK = -tBin^l) ^ |K,r)(K,r^ 



n 



IK 



^ ^(|K,0)(K,r^| + |K,rO(K,0| 
+[/k|K,0)(K,0|, 



where r = IbV^Iti + 1 — \/n{n + 1)], 



— Uib — Ubb 



4, 



iK r' 



(5) 



(6) 



(0,r') 



and 
|K,r> 

|K,0) 



^=^_^e^^-^44^r&i,|0) for r^O 
^JNn{n + 1) R 

^^e'K-%|0)=4|0>. 



(7) 



R 



K is the center-of-mass momentum of the effective two- 
particle system. Since K is a conserved quantity, the 
states with a given K generate an invariant subspace of 
1-L leading to an effective single-particle problem in the 
relative coordinate r. This effective particle moves in a 
K dependent central potential given by V^(r) = /7K(^r,o- 
The hopping amplitude is (n + l)^^ except for bonds con- 
taining the origin where it is equal to ^Jn{n + l)tB- In 
an optical lattice, the impurity atom can be accelerated 
with forces that are invisible to the 6-bosons by using the 
rece ntly demonstrated species-specific dipole potentials. 

One- Dimensional Case. Since the potential V^(r) can 
be attractive, we look for bound-state solutions of Hk- 



(b) 



E = nUr„+Er, 



r = 

FIG. 1: (color online) Lowest energy eigenstates of a MI with 
n bosons per site (n = 3 in the figure) doped with a single 
impurity (black particle), (a) The "particle" and the "hole" 
occupy different positions and the vector r is the relative po- 
sition, (b) The "particle" and the "hole" occupy the same 
position (r = 0). 



For^i: 
as 



1, the exact ground state of can be expressed 



(n + 1) 



= |K,0) + ^e-^/^-|K,r) 



with = n(l - e-^/^^)/(l + n - e"^/^^), and 



,(8) 



.1/Ak _ 



U, 



K 



r/2 



2tB{n+l) ' Y 4t|(n + l)2 + 



The corresponding eigenvalue is 



-2ti3(n+ l)cosh(l/AK). 



(9) 



(10) 



This equation has a physical solution (Ak > 0) for /7k < 
Uc with Uc = -2tB- 



Ak = 



2tBn 



(11) 



Note that the size of the polaron, Aq, diverges for Uq 
Uc. The solutions for [/k > Uc are unbound particle- 
hole states: the impurity- hole item and the boson that 
was promoted to the next Hubbard band propagate in- 
dependently while scattering in each other's vicinity. 

The attractive potential /7k has to reach a critical 
value for the stabilization of the bound state because the 
hopping amplitude is smaller for the bonds that include 
the origin. Since we are assuming that t/,^^ > 0, ejl^ has 
its minimum at K = 0. The effective mass, mj, of the 
exciton is given by the equation: 



1 



l«|K,0)|^ 
mi 



(12) 



where m/ = {2\ti\)~^ is the bare (lattice) mass of 
the impurity in units in which the unit length is given 
by the lattice constant and h = 1. The identity 
([12]) follows from the Hellmann-Feynman theorem and 
Sk('0kI-^' ^)|k=o = 0. By taking the large Aq limit, we 
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obtain the relation between m'' /mi and Aq near the crit- 
ical point: 



mi 



nXo 



(13) 



In the strong coupling regime the effective mass of the 
MI impurity polaron is proportional to its size. This 
behavior, very different from that of lattice polarons in- 
duced by electrons in condensed matter, is caused by the 
unusual mode of transportation: the impurity can hop 
only when the displaced boson and the hole mutually 
annihilate (r=0). As a consequence of this correlated 
motion, the effective mass is proportional to |('^/;q|K, 0)p 
(see Eq.([l2l)). 

To test the validity of H we computed the correlation 
function Cid{'^) between the positions of the impurity 
and the site with one additional boson by solving the 
original Hubbard model with the constraint of no more 
than two bosons per site. We used the Density-Matrix 
Renormalization Group (DMRG) pIBI, fl6| . Figure [2] shows 
the comparison between the numerical results for a chain 
of L = 40 sites and the analytical results given by Eq.®. 
The upper curves correspond to Uib = UBB^tB and the 
lower ones correspond to Uib = UsB — '^tB' As expected, 
both results coincide in the strong-coupling limit. Note 
that the extra boson is displaced over only a few lattice 
sites (r < 6), so that this physics can be realized in to- 
day's optical lattices which typically have a linear size 
corresponding to 100 sites or so. 
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FIG. 2: (color online) Correlation function between the posi- 
tions of the impurity and the site with one additional boson 
for n = 1. The dashed line is the analytical result of Eq.JS]). 
The full lines were obtained by solving 1-L (with a constraint of 
no more than two particles per site) by means of the DMRG 
method in a chain of L = 40 sites. The upper (lower) curves 
correspond to Uib — Ubb + {Uib — Ubb — 2ts). 

General Case. The bound states of 1-L can be found in 
dimension by using the Green's function formalism 
We construct a basis of states that diagonalizes Hik 
and we separate into three terms, = ^ok+'^k + 
Hj^, where 



6±|K,^±)(K,V^^ 



(14) 



e± = (t/K ± Ck)/2, Ck = V/7^ + 4r2, 



IK 



(0,rO 



and u± = (1±/7k/Ck)^^^/a/2- By introducing the Green 
operators 



Ck(^) 



1 



1 



, G+(z) 



1 



(15) 



we obtain the G^(z)-operator from G^{z) by expanding 



in 



G+ = G°c + G°,|K,V+) 



l-e+GO,(^+,^+) 



(K,^+|G°c 



with G^(^+^^+) = (K,V+|G^|K,tA+>. Similarly, by 
considering as the unperturbed Hamiltonian and ex- 
panding in the perturbation H^, we find 



G 



K = G++G+|K,^-)- — i 

i — e Lx^ 



(K,V^-|G+, 
(16) 



with G^(V^",V^") = (K,V^-|G^|K,V^-). The exact dis- 
persion relation of the bound state is obtained from the 
poles of Gk(0,0) = (K,0|Gk|K,0). 

We consider e^, Aq and m* near the critical point, 
Uo = Uc Identifying Ak in d-dimensions from the 
asymptotic behavior of the bound state wave function: 
(V^^l K,r) - ^-r/x^i^{d-i)/2 J. _ (r,0,0,---) and 
r 1, we obtain 

4 = -2tB(n + 1) - 1 + cosh(AK^)] . (17) 

Thus, Ao cx ^^'^ where A5 = —2dtB{n + 1) — eg is the 
composite particle binding energy. From Eq. (p!2]) . we 
obtain the effective mass, m*, 



mi 



(18) 



K=0 



Near the critical point, the large polaron size washes out 
the dependence on specifics that occur on the scale of the 
lattice constant. Hence, we expect the binding energy to 
vanish with the same scaling laws as the binding of a 
single impurity [l^ (2Q| : 



A5(0) (X (/7k - Ucf for d=l and 3, 
C 



A5(0) oc exp 



for d = 2. 



(19) 



where the constant C depends on microscopic details of 
%. These equations lead to the relation between m* and 
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FIG. 3: (a) Kz dependence of spectral weight for the simple 
cubic lattice {Ib — ti — 1^ n — 1^ Kx — Ky — 0, and Ubb — 
UiB —7). Each line is shifted by Kz/n for visualization. We 
use rj — 2~^^ as an infinitesimal constant, (b) Uo dependences 
of the bound state energy Cq (top), linear size of the polaron 
Ao (middle), and ratio of effective and bare masses m* /mi 
(bottom) for n =1. The dashed line indicates Uo — Uc ■ 



Ao: 



— (X Ao for d = 1 and 3, 
mi 

— (X AQ(lnAo)^m/ for d = 2. 

mi 

For the case of a simple cubic lattice in d=3^ Fig. [3] 
shows the spectral density obtained from Eq. (p!6|) . The 
bound states for different values of K that appear be- 
low the bottom of continuum spectrum ([/k < Uc) form 
the polaron band. Uc is obtained by finding a root of 
1 — e~G^(V^~, V^~) = at z = —2dtB{n-\-l) as a function 
of n. Figure ^h) shows Uo dependences of eg, Ao and 
m* /mi for n = 1. Note that the effective mass m* and 



the polaron size increase sharply with decreasing Uc — Uq. 
By superimposing a linear or harmonic impurity-specific 
potential and observing the subsequent acceleration or 
oscillation, cold atom experiments can measure m* di- 
rectly. 

The divergence of the effective mass for Uo Uc is 
cut off at a finite value m^g^^ = Uib / {"^TititB) when the 
next order correction in t^/UiB is included implying that 
m* can be tuned over a large spectrum of values ranging 
from mi = l/2t/ to Uib^ui / {2ntB)' On the other hand, 
the polaron size can be tuned from zero to infinity by 
changing the effective potential Uo or, equivalently, the 
difference between Uib and Ubb- Such dependence of 
m* on the particle-hole pair size is not shared by the 
usual two-particle bound states such as condensed matter 
excitons. In those cases the mass of the bound state is 
well approximated by the sum of the masses. The Mott 
polaron mass increases with size because, to linear order 
in t/, the impurity can only move when the particle and 
hole share the same site (see Figdfb)). This strongly 
correlated mechanism leads to a highly tunable m*. 

The control offered by cold atom technology over po- 
laron size and mass hints at the intriguing prospect of 
studying the quantum phase transition between a "con- 
fined" gas of boson-impurity polarons and two "decon- 
fined" gases of impurities and bosons in the next Hub- 
bard band. The driving parameter of this transition is 
the difference Uib — Ubb- By reducing the ratio UiB/tB^ 
while keeping Uo fixed, it would be possible to vary the 
polaron mass near the transition. The competition be- 
tween potential and kinetic energies could also lead to an 
intermediate crystallization of polarons in the large m* 
regime. 

We are grateful to S. A. Trugman for many useful 
discussions. Work at the LANL was performed under 
the auspices of the U.S. DOE contract No. DE-AC52- 
06NA25396 through the LDRD program. A.E.F. thanks 
NSF for support through Grant No. DMR-0955707. 



[1] M. Greiner and S. Foiling, Nature 453, 736 (2008). 

[2] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[3] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 
[4] M. Greiner, O. Mandel, T. Esslinger, T. Hansch, and 

I. Bloch, Nature 415, 39 (2002). 
[5] M. Greiner, O. Mandel, T. Hansch, and I. Bloch, Nature 

419, 51 (2002). 
[6] G. Campbeh, J. Mun, M. Boyd, P. Medley, A. Leanhardt, 

L. Marcassa, D. Pritchard, and W. Ketterle, Science 313, 

649 (2006). 

[7] W. Bakr, J. Gillen, A. Peng, S. Foiling, and M. Greiner, 

Nature 462, 74 (2009). 
[8] J. Sherson, C. Weitenberg, M. Endres, M. Cheneau, 

I. Bloch, and S. Kuhr, Nature pp. 68-72 (2010). 
[9] W. Bakr, A. Peng, M. Tai, R. Ma, J. Simon, J. Gillen, 



S. Foiling, L. Pollet, and M. Greiner, Science 329, 547 
(2010). 

[10] C.-L. Hung, X. Zhang, N. Gemelke, and C. Chin, Phys. 

Rev. Lett. 104, 160403 (2010). 
[11] S. Inouye, M. Andrews, J. Stenger, H. Miesner, 

D. Stamper-Kurn, and W. Ketterle, Nature 392, 151 

(1998). 

[12] 1. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

[13] J. Catani, G. Barontini, G. Lamporesi, F. Rabatti, 
G. Thalhammer, F. Minardi, S. Stringari, and M. In- 
guscio, Phys. Rev. Lett. 103, 140401 (2009). 

[14] G. Lamporesi, J. Catani, G. Barontini, Y. Nishida, M. In- 
guscio, and F. Minardi, Phys. Rev. Lett. 104, 153202 
(2010). 

[15] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 
[16] S. R. White, Phys. Rev. B 48, 10345 (1993). 



5 



[17] E. Economou, Green's functions in Quantum Physics 

(Springer- Verlag, Berlin, 2006), 3rd ed. 
[18] The word 'novel' here refers not only to the cold atom 

environment but also to the properties of the polaron. 
[19] The increase in polaron size with increasing Uib con- 



trasts with the size of an impurity atom in boson super- 
fluid 

[20] The impurity only affects the diagonal energy of the site 
at the origin. 



